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S . Some recent QCD results from the OPAL Collaboration are summarized. In par- 

^ 1 ticular: a test of color reconnection models and a search for glueballs using gluon 

^~i • jets with a rapidity gap; a study of unbiased gluon jet properties using a new tech- 

vQ ' nique called jet boost algorithm; a measurement of the strong coupling constant 

I as from radiative events. 

K^ ■ 1 Introduction 

mi, 

T^lj- ■ All the following analyses use multihadronic Z" decay events observed with the 

O I OPAL detector ^ at the LEP e+e"conider at CERN. In the first section, results 

>0 . on color reconnection and glueball searches published in a recent paper 3 are 

~r I summarized; section two presents a study j^ of unbiased gluon jet properties and 

^-v . the last section is dedicated to a new measurement [5] of as from radiative events. 

y ■ 2 Gluon jets with a rapidity gap 

Oh! 

O . 2.1 Test of color reconnection models 



^ 



X 



Color reconnection (CR), i.e a rearrangement of the color structure of an event, 
can be a source of rapidity gap events. In the standard Monte Carlo models CR is 
not included, so the color flux tube is stretched from a quark to the corresponding 



^ , antiquark without crossing. If CR is included in the Monte Carlo, we can have 

more complex diagrams, in which string segments can either cross or appear as 
disconnected entities whose endpoints are gluons. In events with an isolated gluonic 
system a rapidity gap can form between the particles coming from the hadronization 
of the isolated segment - often the leading (highest rapidity) part of a gluon jet - 
and the rest of the event. Thus rapidity gaps in gluon jets provide a sensitive means 
to search for color reconnection effects. 

OPAL used this idea to test three CR models: one implemented in the Ariadne 
Monte Carlo 6 , one in Herwig [7] and a model introduced by Rathsman [5| in 
the Pythia Monte Carlo. We refer to these as the Ariadne-CR, Herwig-CR, and 
Rathsman-CR models, respectively. 

All the events are forced to three jets using the Durham 9 jet finder with a 
variable value for the resolution scale jjcut- The two quark (q or q) jets are identified 
using a b-tagging technique [HI] and the remaining jet is the gluon jet. 

The rapidity gap is defined using charged and neutral particles, requiring the 
© Institute of Experimental Physics SAS, Kosice, Slovakia 1 



M. Giunta 



-I— p-i— I— I— r 



-■IZ 



0.015 



0.01 



0.005 



— I— p-i— I— I— p-i— I— i- 
i OPAL 

Jetset7.4 

— ■ Rathsman-CR 
■ ■ • - Ariadne 4.11 

— Ariadne-CR 
Herwig 6.2 
Ti Quark jet 
background 




leading 



Figure 1. Distributions of ri'^^^^^^ in the leading part of gluon jets. "N" represents the total 
number of selected gluon jets and "Ngap" the number of gluon jets with a rapidity gap. 



smallest particle rapidity in a jet ymin > 1-4, or, if this is not the case, the maximum 
difference between the rapidities of adjacent rapidity-ordered particles lS.y„iax > 1.3 
(this second sample is not used for Herwig-CR as explained in 3 ). 

The distributions of charged multiplicity [n'f^^j^^^ ) and total charge [Qieading) 
of the leading part of gluon jets are normalized to the total number of selected jets 
before the rapidity gap requirement. The results for the n'll^adinq distribution are 
shown in Fig. ^ Both Rathsman-CR and Ariadne-CR predict a large excess of 



entries at n 



ch 
leading 



2 and 4 and also at Q 



leadi', 



(figure not shown) . These are 



consequences of events with an isolated gluonic system in the leading part of the 
gluon jet, which is neutral and decays into an even number of charged particles. 
Herwig-CR predicts a less striking excess for 3 < nf^ading — ^ ^^^ Qieading = 0. 

A re-tuning of the Rathsman-CR and the Ariadne-CR models is also tried, but 
in both cases we cannot find a set of parameters allowing the models to describe 



the gluon jets distributions (i.e. 



n 



leading 



and Qieading) whilc Continuing to provide 



a good description of the general features of inclusive Z" decays. It is concluded 
that the Rathsman-CR and the Ariadne-CR models are both disfavored, while it 
is not possible to obtain a definite conclusion concerning the Herwig-CR model. 



2.2 Glueball search 

In |llj . gluon jets with a rapidity gap are also proposed as a potentially favorable 
environment for the production of color singlet bound states of gluons, such as 
glueballs. Therefore, invariant mass spectra in the leading part of the selected 
gluon jets are examined. Since it is a search for anomalous resonant structure, the 
data are compared to the predictions of the models without color reconnection, 
i.e. Jetset, Ariadne and Herwig. These models do not contain glueballs. Three 
different distributions are examined: the total invariant mass of the leading part 
of the gluon jets with Qieading = (since glueballs are electrically neutral), the 
distribution of invariant mass of two oppositely charged particles in the leading 
part of the gluon jets and the corresponding distribution of four charged particles 
with total electric charged zero. No evidence for anomalous production of scalar 
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Figure 2. The ratio between the mean charged particle multipUcities of unbiased gluon and uds 
flavored quark jets, r^/g, as a function of jet energy E* 



particles was observed. 



3 Unbiased gluon jets' properties using the jet boost algorithm 

The Lund theory group proposed a new method, caUed jet boost algorithm |12j 
(BA), to determine properties of unbiased gluon jets. After testing the method, 
OPAL used it to study unbiased gluon jets at different energies and then compared 
the results with theoretical predictions. 

The BA is based on the color dipole model of QCD and can be summarized in the 
following way: in three-jet events, symmetric with respect to the gluon direction, 
the two independent dipoles (one connecting the q and the g and the other going 
from the q and the g) can be boosted and combined to yield the dipole structure 
of a gg event. 

Experimentally, three jets are reconstructed in each multihadronic event and 
the gluon jet is identified with the same technique used in the previous analysis, 
obtaining a final sample of 25396 events with a gluon jet purity of 85%. Then each 
event was boosted to a frame where the angle between the q and the g was the same 
as the angle between the q and the g {9qg — 9qg — 9), yielding a symmetric event. 
The multiplicity of the gluon jet is given by the number of particles lying inside 
the cone-like region defined by the bisectors of 9qg and 9qg and the corresponding 
energy scale of the unbiased gluon jet, E* is equal to the transverse momentum of 
the gluon jet, defined by [Oj : 



PJL. gluon 



(1) 



where Sij {i,j = q, q, g) is the invariant mass squared of the ij pair and s = E'^ ^ . 
The events are divided in 7 subsamples according to the gluon jet energy E*. 

Using Herwig 6.2, we compare results obtained with the BA method on sim- 
ulated three-jet events with results for unbiased gluon jets from simulated color 
singlet gg events. We compare the multiplicity distributions (in the 7 energy in- 
tervals) and we found a good agreement for E* > 5 GeV; we also measured the 
fragmentation functions which agreed for E* > 14 GeV. 
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After establishing the validity of the method, many measurements of unbiased 
gluon jet properties are performed: the mean multiplicity (in the 7 energy bins), 
actually the most precise result for 5.25 < E* < 20 GeV; the first measurement of 
the factorial moments F2 and F^ of unbiased gluon jets multiplicity distributions 
over an energy range; the ratio between the mean charged particle multiplicities of 
unbiased gluon and uds flavored quark jets, Vg/g, as a function of jet energy (see 
Fig. 121). The measurements are compared to many different theoretical predictions 
and an overall global agreement was found. Finally, the fragmentation function was 
measured for jets with E* > 14 GeV and was fitted using the DGLAP evolution 
equation. The fit gives a good description of the measurements and yields a result 
for the strong coupling constant consistent with the world average: 

as{mz) = 0.128 ± 0.008{stat.) ± 0.015{syst.) (2) 

which provides a unique consistency check of QCD. 

4 Measurement of ^5 from radiative events 

Hadronic final states with a high energy isolated photon (e+e"— > qq^) are selected 
and the strong coupling constant as is extracted by fitting event shape variables 
for the reduced center-of-mass energies ranging from 20 to 80 GeV. This procedure 
is allowed by assuming that photons emitted before or immediately after the Z'^ /j 
production do not interfere with QCD processes. 

Our signal is given by multihadronic events with ISR/FSR high energy pho- 
tons {Ej > 10 GeV). The possible backgrounds are: non-radiative multihadronic 
(NRMH) events, two photon processes and r pair production. 

Some isolation conditions are imposed on electromagnetic clusters in order to 
have well isolated high energy photons: the cluster must have an energy larger than 
10 GeV, be located in the barrel region (i.e. | cos{9ec)\ < 0.72) of the detector; the 
angle with respect to the axis of any jet, ajg°, is required to be larger than 25°, the 
sum of tracks momenta and the sum of energy deposition in the electromagnetic 
calorimeter within a cone of 0.2 radiant around the candidate must be both less than 
0.5 GeV. After the isolation cuts 11625 clusters are retained, but the background 
from NRMH events is still large, 52.8%, due to clusters arising from 7r° decay. 

In order to further reduce this background a likelihood ratio method with four 
input variables (| cos{9ec)\, ckjeti duster shape fit variable, distance between elec- 
tromagnetic cluster and associated presampler cluster) is applied. The events are 
divided in seven subsamples according to the cluster energy {Eduster) and the cut 
on the likelihood value is chosen in order to keep 80% of signal events. 

As mentioned in ^1] , JETSET fails to reproduce the observed rate of isolated 
electromagnetic clusters, so the background fraction was estimated directly from 
our data by two independent methods: by fitting the likelihood distributions in 
the data with a linear combination of the Monte Carlo distributions for signal and 
background events or by assuming isospin symmetry (charged hadrons satisfying 
the same isolation cut criteria are selected and the rate of isolated neutral hadrons 
are obtained from the rate of isolated charged hadrons). The final background due 
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to NRMH events is < 10% for ^cluster = 10-35 GeV and 10-15% for ^e/«ster=35-45 
GeV, the two photon processes are < 0.01% and the tt events are 0.5-1.0%. 

After boosting the hadronic system into the center-of-mass frame of the hadrons, 
the following event shape variables are calculated: Thrust (T), Heavy Jet Mass 
{Mh) and Jet Broadening Variables {Bt and Bw)- The normalized background 
distributions are subtracted from data at detector level and then, in order to cor- 
rect for detector effects (acceptance, resolution), bin- by- bin correction factors are 
applied and distributions at the so called hadron level are obtained. The determi- 
nation of as is obtained by fitting perturbative QCD predictions to the event shape 
distributions corrected to the hadron level. The 0(a|) and NLLA calculations are 
combined with the ln(i?) matching scheme and corrected to hadron level multiply- 
ing by a hadronization correction factor obtained using Monte Carlo informations. 
The fit to the event shape variables uses a least x^ method with as{Q) treated as 
a free parameter. The regions with small and uniform corrections are selected as 
fitting ranges. Values of as for each event shape variable and for all event shape 
variables combined are fitted to the solution of the renormalization group equation 
at NNLO described in Jl]. The following value of A]^j„ is derived from the fit: 

AJ^I^^ = 0.2027 ± 0.0141(stsi.)t|]J^^°(sysi.) GeV [x^/d.o.f. = 33.1/6) (3) 

All values of as are then propagated to the energy scale of Mz and a combined 
value of as for all the reduced center-of-mass energy subsamples and all event 
shapes variables is extracted: 

as{Mz) = 0.1176 ± 0.0012(sta<.)tS]:™i(sj/si.) (4) 

This agrees with the previous OPAL analysis using non-radiative LEPl (i.e. 
Ec.m. = 91 GeV) data and the PDG world average. 
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